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Abstract 

It is well known that matter effects affect the way neutrinos oscillate. 
The amplitude of oscillation in matter can be either enhanced, compared 
to the amplitude in vacuum, or suppressed, depending on the density of 
matter at the vicinity of a neutrino. Enhancement is less probable to 
occur than suppression. 

This article demonstrates how matter effects can result into non-oscillating 
neutrinos even in vacuum. 



1 The evolution equation in matter 

We will firstly consider the case of two generations of neutrinos in matter. In 
the basis of flavor, the evolution equation of a neutrino in the 2-Dimensional 
flavor space takes pQ the form: 



d_ ( Ve \ _ f -A + M B 

where A = cos 29 a , B = ^ sin 29 and M = y/2G F N e . N e is the numeri- 
cal density of electrons, Gf is Fermi's constant, Am 2 = rv? v —rr? v and 6g is the 
mixing angle in vacuum, the one and only parameter of the 2x2 mixing matrix 
for 2 generations. At any moment, the probability to observe the neutrino as a 
v e equals |f e | 2 and to observe it as a equals |^| 2 . 

For constant N e , equation JIJ can be solved analytically, as we will see in 
the next paragraph. 



2 Constant N e 

In order to solve eq. JTJ we need to diagonalize the hamiltonian H 



-A + M B 
B A 



It is easy to show that, for constant Ne, the oscillation probabilities are (pQ, 
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where 



and 



P v ,-, Vft = P^„ e = sin 2 26* sin 2 tt— (2) 
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with Ea,b being the eigenvalues of matrix H , which are: 



(4) 



M ± ^/Af 2 + A(A 2 + B 2 - AM) 
Ea.b = g ( 5 ) 

This solution means that neutrinos in matter of constant density oscillate 
in a sinusoidal way, like in vacuum, while the amplitude and wavelength 1 of 
oscillation in matter are different from those in vacuum. The difference results 
from the non-zero M term in equations © and Q. When N e tends to 0, then 
M also tends to and 8 — > 8q and L m — ► L = 4^^^-. 

An important result of this solution is that oscillation amplitude can be 
either enhanced of suppressed with respect to oscillation in vacuum, due to 
matter effects. This may be demonstrated by ranging M in eq. from to oo. 
For M — > +oo we have tan (28) — > =>• 8 — > 90°, which means zero amplitude, 
suppression of oscillation. 

There is only one value of M which maximizes amplitude by giving 8 = 45°. 
This phenomenon is called MSW-resonance 3 and happens when: 

M = Mmsw = 2A^ V2G F N e = cos 28 a (6) 

Regarding the wavelength of oscillation L m , for zero ./V e it equals L$ while 
for N e — * +oo it tends to zero. 



3 Variant N e 

As neutrinos travel through celestial bodies and other objects, they cross regions 
of extremely high and extremely low electron densities. Thus, solving eq. QJ for 
constant N e does not help much. But when N e is a function of position 2 (x), 
then eq. gets extremely complicated and can only be solved numerically 3 . 

x By saying "neutrinos oscillation", what we actually mean is "oscillatory variation of the 
survival/transition probability of neutrinos as time goes by". Thus, amplitude and wavelength 
refer to the characteristics of the sinusoidal function which expresses their survival/transition 
probability. 

2 N e is a function of position, but for neutrinos traveling at (almost) the speed of light, 
it is convenient to replace x with the neutrino time of flight t, using Natural Units where 
c = h = l. 

3 There are only a few exceptions, like the well known adiabatic approximation (I5I.ITIV 
where the evolution equation with variant electrons density can be given an analytical ap- 
proximate solution 
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Figure 1: (a): The density profile for which we solved the evolution equation. 
It is a gaussian function which most of the time exceeds Nmsw- (b): The 
resulting survival probability of a neutrino which initially was a pure v e . We 
can see the region of suppression at the center of the density distribution. We 
have assumed E=3 GeV, 9 = 32°, Am 2 = 7.2 • 1(T 5 eV 2 . 

Either numerically or analytically, solutions of eq. Q give us the amplitudes 
v e {t) and Vn(t) as functions of time. Those amplitudes are complex numbers, 
but we usually express the solution in terms of |i/ e (i)| = P {y e \ t) and |f M (i)| 2 = 
P [v^i\ t), as done in eq. J2J), because we have a more intuitive understanding of 
those probabilities. However, as it will have been explained by the end of this 
article, the complex nature of these amplitudes is important and my lead to 
observable phenomena. 

When N e varies, P {v e ;t) and P(v^\t) may vary in a chaotic fashion, de- 
pending on the density profile N e (t), the initial conditions v e ,ij,(t = 0) 5 the 
energy E of the neutrino and the input parameters Oq and Am 2 . An example of 
such an irregular oscillation can be seen in Fig.^ Nevertheless, in space (time) 
intervals where iV e does not change very rapidly, the approximation of constant 
N e holds well and oscillation patterns in those regions look like those of eq. J5J . 
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Another thing that should be noticed is that MSW-resonance may occur 
instantly if the density profile has not been adjusted so as to satisfy eq. © for 
long intervals. On the other hand, suppression is observed for N e ^> Nmsw, 
which holds for a whole range of density values, not just for a specific N e . 
So, if things have not been intentionally set up to be different, suppression is 
a phenomenon which may last for longer and is more probable than MSW- 
resonance. 



4 The complex nature of v e ^ 

Equation (JTJ is a system of complex differential equations and its solutions are 
complex functions of time, as stated earlier. The real and imaginary parts of 
u e {t) and Vfi{t) affect the real and imaginary parts of v e (t + dt) and v^it + dt). 

As for every differential equation, the solution of eq. |QJ depends on the initial 
condition v e ,n(t = 0), which involves v e (Q) and ^(0) as complex numbers: 

v a = \v a \ e l0Q , a = e,/x (7) 

Starting from the same \u a \ and different <p a leads to completely different solu- 
tions P(v a ,t), even in vacuum 4 . However, people don't pay much attention to 
the arguments <j> a . The reason is that neutrinos are always produced in weak 
vertices. All experiments show that weak bosons (W^^Z ) couple with neu- 
trinos in eigenstates of flavor. Thus, all neutrinos start their journey from the 
initial condition: 

P(v e ;0)\ ( 1\ _ f0' 



P(^;0)J \0J \1 
equivalently: 



I V o 



or m (8) 



Of course, initial condition (|SJ| says nothing about the arguments <j) a . They 
can take any value. But if they are arbitrary and also drastically affect P(y a \ t), 
then experiments should have observed extremely messy oscillations, different 
for identical neutrinos with only different <j) a . The reason that this is not what 
happens is that the impact of <p a on P(u a ;t) is canceled when eq. JHJ holds, 
which means always. It can be proved analytically and also be tested numeri- 
cally. 

In brief, neutrinos are eigenstates of flavor, not of mass, and that is why 
they oscillate. That is also why <fi a have no observable effect. 



5 The condition of non-oscillation in vacuum 

There are at least two methods to find the condition of non-oscillation of neu- 
trinos in vacuum. The brute one is to solve analytically eq. Q with the most 

4 As we will see soon, there is only one case where <j> a makes no difference, and this case is 
the one nature always prefers. 
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general expression of initial conditions, find the expression for P(y a \ t) in terms 
of the initial conditions and then demand that ap ^'°'^ = Vt. This demand 
imposes a condition to the initial conditions which, if satisfied, guarantees con- 
stancy 5 of P{v a ; t). 

Constancy of P{y a ;t) is also guaranteed if the neutrino initially is in an 
eigcnstate of mass, because E = y^p 2 + m 2 so, assuming definite momentum 
p, definite mass means definite energy E, which means constant state. So, 
the second, equivalent but much more elegant method to find the condition of 
non-oscillation is to demand that initially the neutrino has definite mass. 

For two generations, the mixing of flavor and mass eigenstates is given by 
the expression: 

v e \ _ / cos 9q sin Oq \ / v\ 
Vu) ~ I - sin 9 cos 8 J U 2 



v\ \ I cos uq — sin vq \ i v e 
v 2 J \ sin 6»o cos 6 I \ v„ 



0) 



with 



v e \ 2 + |^| 2 = 1 (10) 



Demanding that initially vi =0, 11U|) and I© give the condition 



tan 6*o and \v e 



\/l-|-cot 2 8 



Demanding that initially v 2 — 0, from i|l(J|) and © we find: 



— = — cot Qn and 



\/l+tan 2 8o 



(11) 



(12) 



Conditions (|11|) and l|12|) do not only refer to the amplitudes of v e ^ , but also 
to the relationship of their arguments <p a , so the complex nature of v e and 
can not be marginated and consider only \v a \. For (|TT1) to be satisfied, v e and 

must be in the same direction on the complex plane (homoparallel), while 
l|12|) demands them to point in opposite directions (antiparallel) . 8q is the only 
parameter in both non-oscillation conditions. Fig.[21shows the configurations of 
u e and on the complex plane, which satisfy the non-oscillation conditions. 

Let's examine how v e and behave as complex numbers, when in vacuum. 
Let me substitute v e and with the symbols x and y respectively. Analytical 
solving of eq. QJ for M = gives x and y revolving on ellipses on the complex 
plane, like in Fig. |3J Those ellipses always have their centers at the point 
where axes intersect (point 0) and are inclined with respect to the axes. Their 
inclinations depend on both xq and ya. The analytical expression of those 
ellipses has been calculated to be: 



Re [x(t)\ = ynf+TTfcos (st - arctan (nfj) 
Im [x(t)} = yiTfTTTf cos yst + arctan ( §^ 



(13) 



5 Constancy of P(y a \t) does not mean constancy of the solution u a (t), but only of |^ Q (t)|. 
In the complex plane, u e (t) and Vu{t) always rotate. 
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Figure 2: Complex pairs like (xo,yo) and (xo,yo), pointing in any direction, 
satisfy (|1 l|l and 1)12(1 respectively, assuming Oq = 32° which determines the radii 
of the circles, xq and yo stand for v e and v^. 



Re [y(t)} = y/E* 
Im [y(t)} = ^Bl 



where 



B\ cos ( st 



(st — arctan fjjf'JJ 



B\ cos ( st - 



arctan 



(14) 



III = XnRe [Ci] + X 21 Re [C 2 ] , B 1 = X 12 Re [Ci] + X 22 Re [C 2 ] 
n 2 = X 21 Im [C 2 ] - X n Im [d] , B 2 = X 22 Im [C 2 ] - X 12 Im [d] 
H 3 = Xulm [Ci] + X 21 Im [C 2 ] , B 3 = X 12 Im [d] + X 22 Im [C 2 ] 

n 4 = X n Re [Ci] - X 21 Re [C 2 ] , B 4 = X 12 Re [Ci] - X 22 Re [C 2 ] 



(15) 



where 



where 



and 




d 



C< _ Xg-XnCl 

° 2 - 



(i+(^)T 1/2 (^ 



2\-V2 / 1 



4£ 



(16) 



(17) 



(18) 



G 



From l|13l) and (|14f) it is obvious that x(t) and y(t) revolve with the same 
angular frequency s, so, if they are not initially collinear in the complex plane, 
they will never be. On the other hand, if we initially select xq and yo to satisfy 
one of the non-oscillation conditions, then the ellipses take the shape of the black 
circles of Fig. 3 and x and y keep revolving coherently, being always collinear. 
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Figure 3: Trajectories of x(t) and y(t) in the complex plane when in vacuum, 
with initial conditions xo — 0.5e l ^,z/o = 0.866e l0 (see the dots on the ellipses) 
and 9o = 32°. Red ellipsis is the trajectory of x(t) and blue is of y(t). Black 
circles are the same we had in Fig. [21 



6 Matter effects in complex trajectories 

Non-oscillation in vacuum would only be fiction, if matter effects were not capa- 
ble of achieving the non-oscillation condition, because nature produces neutrinos 
obeying eq. ijHJl, which contradicts and (|12J) . 

Numerical solution of eq. Q with variant N e has revealed that the trajec- 
tories of x and y in the complex plane get extremely complicated (see Fig. 
and can happen to satisfy one of the non-oscillation conditions, even instanta- 
neously. In other words, a neutrino in matter may transiently get into a mass 
eigenstate. 

From the left part of Fig. Q]we realize that, as density increases and then 
decreases, x(t) and y(t) drift away from the ellipsoidal trajectories they initially 
followed in vacuum. When density becomes zero again, they return to elliptical 
trajectories which are different from those they initially followed when in vac- 
uum. The reason is that the neutrino returned to vacuum from different initial 
conditions, i.e. the x and y it had at the moment it exited the matter region. 
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Figure 4: Right: The survival probability of a neutrino which initially was a 
v e with xq — 1 • e l s , yo = and crosses a region of high matter density, such 
as a star. We have assumed a gaussian N e (t) of the form of eq. fTT?)) with p = 
10- 17 GeV 3 , t = 2 ■ 10 24 GeV- 1 and cr = 10 23 GeV- 1 . Suppression is obvious 
around t = t , as N e highly exceeds Nmsw- The oscillation is drastically 
modified by the passing through the star. Left: The complex trajectories of 
x(t) (red) and y(t) (blue). Before entering the star, x(t) circulates along the red 
ellipsis (red and blue dots indicate xq and j/o respectively) and y(t) along the 
blue ellipsis which is degenerated into a line for those (naturally default) initial 
conditions. 

The most critical interval is while the neutrino lies in the middle of the matter 
region, where its oscillation is highly suppressed. Then, y(t) moves along the 
curly blue line, around y(t) ~ something + i -0.9. At the same time, x(t) rapidly 
spins along the red circles of radius ~ 0.3. In what trajectories they will end 
up depends on when the spinning of x(t) and the curling of y(t) will stop, i.e. 
when N e will decline. 

If iV e rises appropriately and diminishes at the right time, then it is possible 
to make x(t) and y(t) satisfy (or almost satisfy) one of the non-oscillation con- 
ditions. After methodically trying several initial conditions (always obeying eq. 
JHJ) and several N e shapes, I found one example of such a 'flavor lens' 6 . The 
neutrino was assumed to initially have xo = 1 • e l ° , yo — and energy E=3 GeV. 
It has also been assumed that 9q = 32° and Am 2 = 7.2 x 10~ 5 eV 2 . 

The density profile found can be seen in Fig. 5, together with the oscillation 
of survival probability of a v e passing through it. N e (t) rises according to the 

e Actually, trial-and-error method would have very hardly lead to any findings. I had to 
implement a 'trick': At each time step of the numerical solution I would check if any non- 
oscillation condition were almost satisfied. If it were, then I would 'cut' the density distribution 
and set it equal to zero from there on. In this way, the neutrino would be dropped back to 
vacuum on the right time. 
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1x10" 2x10" 3x10" 4x10" 5x10" 

Figure 5: A density profile which acts as a flavor lens for v e of E=3 GeV, 
provided that O = 32° and Am 2 = 7.2 • 1(T 5 eV 2 . 

function: 2 

JV e (f)=A)exp(- (f ~J, o) j (19) 

with 

to = 4 - 10 24 GeV- 1 
po = 1CT 16 GeV 3 
a = 10 24 GeV' 1 

N e (t) = for t > 3.03 • 10 24 GeV- 1 

The fact that I found this specific density profile does not mean that it is 
the only 'flavor lens' for this kind on neutrinos. The same distribution, with 
its center to translated by an integer number of vacuum wavelengths Lq, would 
have the same result. Whether a density profile is a flavor lens or not depends 
on its position with respect to neutrino's production point, on neutrino's energy, 
on its shape and on the initial conditions of the neutrino (xq, yo)- 

It has been tested and confirmed that 4> a do not affect the oscillation through 
the flavor lens, as long as JSJ was initially true. In the same way that (JSJ 'hides' 
the impact of <fi a in oscillations in vacuum, it also also prevents the role of <fr a 
from being revealed when matter is present. So, the same oscillation would 
occur whether (xo = 1 • e M , yo = 0) or (xq = 1 ■ e lTr / 3 , yo = 0) etc. (see Fig. 8). 

7 3 Generations 

In 3 generations things do not differ much in principle. The equation of evolution 
is a system of 3 differential equations, where there is a matter term equivalent 
to M we had in 2 generations [Q. 

Demanding that a neutrino is at an eigenstate of mass i/£, instead of an 
eigenstate of flavor i> a , we derive the non-oscillation condition in 3 generations. 
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Figure 6: Survival probability for a v e of E=3 GeV, as it passes through the 
flavor lens of Fig. 03 When it has emerged from the matter region it is almost 
non-oscillating. The fact that it undergoes very slight oscillation is because the 
non-oscillation condition was not exactly satisfied when matter dropped to zero, 
but was almost satisfied. 
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Figure 7: The trajectories of x(t) (red) and y(t) (blue) as a 3 GeV neutrino 
with xq = 1 • e l °, yo = passes through the flavor lens of fig. El 
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Figure 8: The same like in fig. but with initial conditions xq — 1 • e l ~ , yo = 0. 
Both trajectories are rotated by tt/3, but their behavior is the same as fig. 
so the survival probability remains the one in fig. |H1 



Starting with the notation: 

3 



v a = ^2 U <xiVi =>Vi= ^2 U* al v a (20) 

i— 1 cx—e,fi,r 

and respecting the completion condition: 

We\ 2 + + Wr\ 2 = 1 (21) 

we find that to have an eigenstate of mass, i.e. to have Vyt = 1, ity = Vi = 0, the 
condition is: 

= l/^/lT^I 2 + \B? 

Av" (22) 



where 





' K 


] 












u*, L 

er 



r; " B=(-2kA-?£) ,2,-n 



In 3 generations, the non-oscillation condition is triple, one for each eigen- 
state of mass, and depends only on the tree mixing angles of the PMNS matrix 
U. 

Non oscillation condition in 3 generations is much more demanding than 
in 2 generations. It demands three complex numbers to be collinear in the 
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complex plane, to point in the right directions and to be of the right sizes. 
However, achieving its satisfaction with an appropriate density profile must not 
be impossible, though it is not easy to find an example. 



Conclusions 

Matter effects in neutrinos oscillation have been discussed. The complex char- 
acter of f M and v e has been emphasized, as density N e affects their trajectories 
in the complex plane, potentially leading neutrinos, transiently, to mass eigen- 
states. 

It has been demonstrated how an appropriate N e (t) can transform a regular 
v e into one which does not oscillate in vacuum. 

Finally, the non-oscillation condition in 3 generations has been presented. 
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